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SUMMARY 

A  mathematical  representation  of  vertical  acceleration  response  characteristics  of  a  helicopter  in  hover 
has  been  developed,  including  blade  flapping,  inflow,  and  rotor  speed  dynamics.  A  Maximum  Likelihood 
parameter  estimation  technique  has  been  applied  to  assess  the  adequacy  of  the  model,  and  to  identify  the 
relevant  parameters,  using  flight  data  from  a  Sea  King  Mk  50  helicopter.  A  number  of  conclusions 
related  to  the  validity  of  the  modelling  approach  have  resulted  from  comparisons  between  predicted  and 
identified  parameters,  and  further  investigation  of  some  aspects  is  indicated. 
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NOTATION 


A  System  matrix 

B  System  matrix 

Cj  Thrust  coefficient 

e  Flapping  hinge  offset  distance 

g  Acceleration  due  to  gravity 

Ib  Rotor  blade  moment  of  inertia  about  hinge 
lz  Total  moment  of  inertia  of  main  rotor  about  the  shaft 

J(£)  Cost  functional 

Ki  Pitch  -  flap  coupling  factor 

m  Helicopter  mass 

mb  Mass  of  a  main  rotor  blade 

Ma  Total  aerodynamic  moment  about  hinge 

Mb  Rotor  blade  moment  of  mass  about  hinge 

Mu  Apparent  additional  air  mass  factor 

n  Measurement  noise  vector 

N  Number  of  blades  of  the  main  rotor 

P  System  matrix 

Total  aerodynamic  torque 
Shaft  torque  of  engines 
R  Covariance  of  residuals  matrix 

R  Radius  of  main  rotor  blade 

t;  ilh  discrete  time  point 

Ta  Total  aerodynamic  thrust 

u  Control  input  vector 

w  Vertical  velocity 

x  State  vector 

z  Measurement  vector 

zj:  Output  vector  estimate 

Z  Vertical  acceleration 
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NOTATION  (cont) 


Blade  flapping  angle 
Blade  Lock  number 
Perturbation  in  thrust  coefficient 
Sensitivity  matrix 

Perturbation  in  non  dimensional  inflow  velocity 
Ratio  of  flapping  hinge  offset  to  rotor  radius 
Rotor  blade  pitch  angle 
Inflow  velocity 

Non  dimensional  inflow  velocity 
Reference  non  dimensional  inflow  velocity 
Unknown  system  parameter  vector 
Rotor  angular  velocity 


1 .  INTRODUCTION 


Because  of  the  complex  nature  of  rotorcraft  aerodynamics  and  flight  dynamics,  any 
global  model  which  adequately  represents  dynamic  flight  characteristics  will  contain  a 
large  number  of  parameters  based  on  numerous  simplifying  assumptions.  Validation  of 
such  models  against  flight  measurements  is  a  difficult  and  time  consuming  exercise.  One 
such  attempt  [1]  with  the  Sea  King  Mk  50  helicopter,  demonstrated  that  while  good 
overall  representation  can  be  achieved,  significant  deficiencies  remain  in  specific  areas.  A 
better  understanding  of  these  deficiencies  can  be  gained  by  isolating  the  particular  aspect 
of  concern  and  focussing  on  the  development  of  an  adequate  model  structure  which  can 
be  more  readily  verified.  The  simplified,  special  purpose  models  which  result,  are  also 
especially  useful  in  studies  related  to  flying  qualities,  stability  augmentation  system 
design,  and  validation  of  in-flight  and  ground  simulators  [2], 

Parameter  estimation  techniques  developed  over  the  last  two  decades  have  been  used 
widely  by  the  fixed  wing  community,  and  more  recently  are  being  applied  in  studies  of 
rotorcraft  flight  dynamics  [3-6].  For  the  latter,  as  well  as  the  problems  of  a  noisier 
measurement  environment,  particular  difficulties  are  presented  by  the  lack  of  a  well 
defined  model  structure,  and  the  presence  of  a  large  number  of  unknown  parameters. 
Nevertheless,  parameter  estimation  techniques  provide  a  useful  tool  for  investigating 
which  features  are  more  important  to  retain  for  an  adequate  model,  and  to  validate 
theoretical  estimates  of  the  corresponding  parameters.  The  probability  of  successful 
application  of  these  techniques  can  be  enhanced  by  making  full  use  of  any  a  priori 
information  available. 

Vertical  axis  response  characteristics  are  important  in  defining  flying  qualities, 
especially  in  flight  close  to  the  ground  involving  bob-up  manoeuvres  [7].  A  simple 
modelling  approach  as  used  in  Reference  [1]  is  not  adequate  for  predicting  transient 
response  of  collective  inputs,  including  the  initial  peak  response  and  subsequent  periodic 
behavior.  In  particular,  it  has  been  noted  [7]  that  blade  flapping  dynamics  can  have  a 
significant  influence  on  helicopter  flight  dynamics  and  handling  qualities.  Similarly, 
Reference  [8]  has  shown  that  dynamic  inflow  can  play  a  key  role  in  the  initial  overshoot 
of  the  vertical  acceleration  to  an  abrupt  collective  input.  Finally,  the  importance  of  engine- 
rotor  dynamics,  and  hence  rotor  speed  changes,  on  vertical  axis  handling  qualities  has 
been  emphasized  in  Reference  [9],  This  suggests  the  desirability  of  including  rotor  speed 
dynamics  in  any  model  of  vertical  axis  response  to  collective  inputs. 

This  paper  develops  a  model  for  vertical  response  in  hover  incorporating  each  of  the 
above  features,  and  by  use  of  time  domain  parameter  estimation  techniques,  in 
conjunction  with  available  flight  data  for  the  Sea  King  Mk  50  helicopter,  establishes  the 
need  to  retain  them  in  an  adequate  model.  Secondly,  parameter  values  extracted  from 
measured  responses  to  collective  step  and  pulse  inputs  are  compared  with  theoretical 
estimates.  Practical  difficulties  due  to  limited  data  are  discussed,  and  some  suggestions 
made  for  further  model  refinement  and  validation. 
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2.  MATHEMATICAL  REPRESENTATION 


Equations  for  vertical  motion,  blade  flapping,  inflow  and  rotor  speed  dynamics  were 
derived  using  the  following  assumptions: 

(i)  Rigid  rotor  blade  with  linear  twist  and  uniform  section 

(ii)  Pure  vertical  motion,  i.e.  all  angular  rates,  fore/aft  and  sideways 
velocities  are  zero 

(iii)  Uniform  inflow 

(iv)  Cyclic  flapping  coefficients  are  zero 

(v)  Incompressible,  attached  flow  over  blades  with  tip  loss  factor  of  1 

(vi)  Changes  in  fuselage  download  neglected 

(vii)  Flapping  and  inflow  angles  small 

(viii)  Only  first  order  terms  in  hinge  offset  retained 

2.1  Vertical  Motion 

The  equation  for  vertical  velocity,  w  (positive  down),  can  be  written  as  follows: 

mw  -  mg  -  Nl^bR  (l-e)P  =  -TA  (1) 


where  m  is  the  helicopter  mass,  mb  is  the  blade  mass  and  R  the  rotor  radius,  N  is  the  total 
number  of  main  rotor  blades,  e  is  the  ratio  of  the  hinge  offset  to  rotor  radius,  e/R,  P  is  the 
flapping  angle  (positive  up),  and  Ta  is  the  total  aerodynamic  thrust. 

2.2  Blade  Flapping 

The  equation  for  blade  flapping  dynamics  is: 

IbP  +  £22(Ib  +  eMg)P  -  Mbw  =  Ma  (2) 


with  Ir  and  Mb  the  blade  moment  of  inertia  and  moment  of  mass,  respectively,  about  the 
flapping  hinge,  Q  is  the  rotor  angular  velocity  and  Ma  is  the  aerodynamic  moment  about 
the  flapping  hinge. 

2.3  Inflow 

The  inflow  dynamics  equation  follows  the  derivation 

Mil  ^5.+  2 ffi  -  JSL  +  l  A 

n  dt  [  qr  3  n 

where  u  is  non  dimensional  inflow  velocity,  Cj  is  the  thrust  coefficient  and  M  i  j  is  related 
to  the  apparent  additional  air  mass  introduced  by  Carpenter  and  Fridovich  [10j. 


of  Chen  and  Hindson  [8] 
Wor  (3) 


Values  of  Mu  vary  from  0.849  in  Reference  [10]  to  a  'corrected'  value  of  128/7571 
(0.543)  in  the  dynamic  inflow  theory  of  Pitt  and  Peters  [11].  The  latter  also  imply  a 
simpler,  linear  alternative  to  Equation  (3) 

Mn  dAu 

— - +  2  \)Rcf  At)  =  ACt  (3a) 

where  Au  and  ACt  are  perturbations  in  inflow  and  thrust  coefficient. 

2 . 4  Rotor  Speed 

The  equation  for  the  rotor  speed  variations  can  be  written 

Izn-2N(lB  +  eMB)(Jpa  =  QH  +  QA  (4) 


where  Qe  is  the  shaft  torque  provided  by  the  engines,  Qa  is  the  total  aerodynamic  torque 
about  the  shaft,  lz  is  the  total  rotor  moment  of  inertia  about  the  shaft. 

The  force  and  moment  terms  on  the  right  hand  side  of  the  above  equations  were 
derived  using  a  simple  strip  theory  detailed  in  Reference  [12].  In  order  to  facilitate 
comparisons  between  model  and  flight  data  via  parameter  estimation,  the  equations  (1) 
to  (4)  are  linearized  about  a  reference  state  of  steady  hover  with  zero  vertical  velocity. 
The  linearized  small  perturbation  model  has  the  form 
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(5) 


The  input  vector  in  (5)  includes  changes  in  collective  pitch,  9c,  and  engine  torque. 
A  priori  values  for  each  of  the  matrix  elements  can  be  obtained  from  the  expressions 
summarized  in  Appendix  A.  Relevant  physical  data  used  for  the  estimating  a  priori  values 
for  the  Sea  King  Mk  50  helicopter  is  given  in  Table  1 
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Table  1  -  Sea  King  Physical  Characteristics 


Rotor  radius  (ft),  R 

31 

Blade  chord  (ft),  c 

1.52 

Hinge  offset  (ft),  e 

1.05 

Blade  twist  (deg),  0t 

-8 

Pitch  -  flap  coupling,  Ki 

0.08 

Blade  mass  (lb),  mb 

181 

Angular  velocity  (rad/s),  Ci 

21.89 

Lock  number,  y 

11.51 

I7 

Ratio  of  rotor/blade  inertia,  j~ 

6 

3.  FLIGHT  DATA 

A  flight  test  program  using  the  Sea  King  Mk  50  helicopter  has  been  described  in  detail 
in  Reference  [13],  and  has  provided  an  extensive  data  base  for  both  performance  and 
flight  dynamics  characteristics.  Data  was  recorded  in  12  bit  form  at  a  sampling  rate  of 
60Hz  and  subject  to  a  range  of  post-processing  procedures  to  reduce  random  noise  levels, 
correct  various  error  sources  such  as  'drop-outs'  and  time  delays,  and  to  ensure  kinematic 
consistency.  The  data  of  particular  relevance  here  are  the  vertical  acceleration  responses  to 
collective  step  and  pulse  inputs  in  hover.  The  cases  selected  for  the  current  investigation 
are  summarized  in  Table  2,  including  information  on  input  type  and  all  up  weight 
(AUW),  which  affects  directly  or  indirectly  several  of  the  parameters  listed  in 
Appendix  A.  Unfortunately  other  input  shapes,  such  as  doublets,  3211,  or  sweeps  were 
not  available.  Nevertheless,  the  table  shows  that  there  were  nine  cases,  each  between  five 
and  eight  seconds  duration,  included  in  the  analysis.  The  first  six  cases  being  at  a  higher 
AUW  of  around  18500  lb  and  the  remainder  at  a  lower  AUW  of  about  16500  lb.  In  all 
cases  the  altitude  was  close  to  sea  level.  For  each  case,  measurement  records  were 
available  for  collective  input  and  vertical  acceleration  response,  as  well  as  engine  torque 
and  rotor  speed,  but  no  information  was  available  for  the  inflow  or  blade  flapping  angle 
variations. 
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Table  2  -  Flight  Cases 


Case 

Collective 

Input  Type 

AUW 

(lb) 

1 

Step  up 

18925 

2 

Step  up 

18835 

3 

Step  down 

18719 

4 

Step  up 

18487 

5 

Step  up 

18449 

6 

Step  down 

18335 

7 

Step  down 

16535 

8 

Step  up 

16590 

9 

Pulse 

16500 

4 .  PARAMETER  ESTIMATION  METHOD 

The  linear  system  given  by  Equation  (5)  has  the  general  form 

Px  =  Ax  +  Bti  (6) 

where  the  vectors  x  and  u  are  the  state  and  input  respectively,  while  matrices  P,  A,  and  B 
contain  the  system  parameters.  The  output  at  time  tj  is  modeled  as 

z(tj)  =  C  x(q)  +  D  u(tj)  +  n(tj)  (7) 


where  n(tj)  is  the  measurement  noise  vector  and  C  and  D  are  also  system  matrices. 

In  the  present  formulation,  the  state  equation  (6)  is  assumed  free  of  process  noise 
while  the  measurement  noise  is  assumed  to  be  zero  mean  white  Gaussian  noise.  The 
maximum  likelihood  (ML)  of  the  unknown  system  parameters,  is  that  value  of  ^  that 
maximizes  a  likelihood  function,  defined  as  the  probability  of  obtaining  the  outputs.  /., 
given  the  parameters,  1;.  Details  of  the  method  can  be  found  in  standard 
references  ( 14,15],  and  amounts  to  the  minimization  of  a  cost  functional 

1  N  l 

J(5)  =  7.Z  [*<«i)  -  ytO]7  R  [*(t0  -  K^(ti)  ]  (8) 
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with  an  estimate  of  R  being  given  by 


1  N 

R=jq-  I  [  '-(>.)  -  z^Ui)]  [/-(ti)  -  ]T 


(9) 


where  z^(tj)  is  the  estimate  of  the  output  at  time  tj.  A  modified  Ncwton-Raphson 
algorithm  is  used  to  achieve  the  minimization  iteratively.  A  measure  of  the  accuracy  of  the 
estimates  is  given  by  the  Cramer-Rao  bound  [  14]: 

Covariance(^)  = 


l[V4(/.^(ti»]TR''  V^tj)) 


t=i 


(10) 


where  the  gradient,  V,  term  is  the  sensitivity  matrix  obtained  from  the  system  equations. 
The  Cramer-Rao  bound  is  a  lower  bound,  and  Reference  |15]  points  out  that  it  is 
significantly  smaller  than  the  scatter  of  the  results.  Nevertheless,  it  is  a  reasonable 
indicator  of  relative  accuracy  of  various  parameters  from  different  cases.  Correlations 
between  parameters  can  also  be  derived  from  the  covariance  matrix  (10) 

A  particular  advantage  of  the  output  error  method  is  that  not  all  the  outputs  need  to  be 
measured  for  parameter  estimates  to  be  obtained.  For  the  present,  only  the  vertical 
acceleration  and  rotor  speed  are  available  for  matching,  while  the  total  number  of 
parameters  indicated  in  Equation  (5)  is  28.  With  the  limited  information  available  and  the 
non-optimal  input  shapes,  it  is  not  possible  to  obtain  reliable  estimates  of  all  28 
parameters,  many  of  which  are  highly  correlated,  even  if  convergence  of  the  algorithm 
were  achieved.  In  order  to  reduce  the  total  number  of  unknown  parameters,  the  a  priori 
information  has  been  used  in  two  ways.  Firstly,  an  order  of  magnitude  analysis  was 
undertaken  to  establish  which  parameters  were  likely  to  be  the  most  significant.  As  a 
result,  several  of  the  less  important  parameters  were  fixed  at  their  a  priori  values. 
Secondly,  many  of  the  remaining  parameters  were  closely  related,  differing  from  one 
another  only  by  a  constant  multiple.  By  maintaining  these  constraints  during  the  parameter 
estimation  process,  the  total  number  of  independent  variables  can  be  reduced 
considerably.  In  Appendix  A,  with  fixed  and  constrained  parameters  indicated  in  the  tight 
hand  column,  the  remaining  number  of  independent  parameters  has  been  reduced  to  14, 
of  which  five  relate  to  the  torque  equation,  and  three  relate  to  each  of  the  other  equations. 
A  priori  values  were  used  as  a  guide  to  initial  estimates  for  the  parameters  during  the 
identification  process.  Even  so,  convergence  was  marginal  when  all  14  parameters  were 
simultaneously  identified,  and  it  was  necessary  to  vary  the  starting  values  of  several  of  the 
parameters,  particularly  in  the  rotor  speed  equation,  in  order  to  achieve  convergence. 
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Because  of  the  convergence  difficulties,  an  alternative,  more  robust  identification  strategy 
was  developed  by  decoupling  the  rotor  speed  equation  from  the  rest  of  the  system 
equations,  as  follows: 

1 .  Set  to  zero  all  parameters  in  the  rotor  speed  equation  except  for  £2{j,  fig,  and 
Qq.  Match  the  rotor  speed  and  vertical  acceleration  records  for  initial 
estimates  of  £1,  z,  p,  and  t)  parameters.  Initial  fits  are  quite  good. 

2.  With  the  z,  P,  and  o  parameters  fixed  at  the  values  obtained  at  step  1,  obtain 
an  improved  £2  match  with  additional  parameters  £ip,  Qw,  and  all  constrained 
parameters,  included. 

3.  Use  the  £1  record  from  step  2  as  an  input  to  the  system,  to  obtain  an 
improved  vertical  acceleration  march  and  revised  z,  p,  and  l>  parameters. 

Steps  2  and  3  above  could  be  iterated  further  if  necessary  but  in  practice  this  does  not 
produce  any  significant  improvement,  either  in  terms  of  fit  errors  or  parameter  accuracy  as 
measured  by  the  Cramer  -  Rao  bound.  The  three  step  procedure  was  successful  in  all 
cases  tried,  and  the  results  were  in  good  agreement  with  those  obtained  with  all  14 
parameters  simultaneously  identified,  when  convergence  was  achieved  for  the  latter. 

5.  RESULTS 

S.l  Time  Histories 

The  effect  of  including  flapping  and/or  inflow  dynamics  in  the  vertical  response  model 
is  illustrated  in  Figure  1,  which  shows  the  best  match  for  vertical  acceleration  achievable 
by  each  model  against  the  flight  data  of  Case  1.  In  Figure  1(a)  the  model  includes  only 
inflow  dynamics  while  Figure  1(b)  has  flapping  but  not  inflow  dynamics  represented. 
Figure  1(c)  has  both  flapping  and  inflow  dynamics,  but  in  none  of  the  cases  has  any 
account  been  taken  of  rotor  speed  variations. 

While  Figure  1(a)  illustrates  that  inflow  dynamics  cannot  reproduce  the  initial 
acceleration  peak,  it  is  clear  from  Figure  1(b)  that  flapping  dynamics  on  its  own  is  also 
inadequate.  However,  with  both  effects  included  an  excellent  match  with  the  initial 
acceleration  response  to  a  step  input  can  be  achieved,  as  seen  in  Figure  1(c). 
Nevertheless,  the  longer  term  periodic  variations  cannot  be  reproduced  without  taking  into 
account  rotor  speed  dynamics.  When  the  complete  model,  as  described  in  Section  2,  is 
used,  an  excellent  match  of  both  vertical  acceleration  and  rotor  speed  records  can  be 
achieved.  This  is  illustrated  in  Figures  2  to  4  for  three  different  cases.  Figure  2 
(Case  4)  is  the  response  to  a  collective  step  up  (note  that  positive  acceleration  is 
downwards).  Figure  3  (Case  7)  is  for  a  collective  step  down  input,  while  Figure  4 
(Case  9)  illustrates  the  response  to  a  collective  'pulse'  type  of  input.  For  the  step 
response  cases  it  is  clearly  possible  to  distinguish  the  initial,  high  frequency,  response 
from  the  longer  term  dynamics,  but  the  time  history  for  the  pulse  input  is  a  combination  of 
both  high  and  low  frequency  dynamics.  While  band  pass  filtering  could  be  used  to 
separate  these,  this  has  not  been  found  necessary  in  this  investigation,  since  the 
uncoupling  procedure  outlined  in  the  previous  section  largely  achieves  the  same  end.  The 
clear  separation  between  short  and  long  term  effects  is  probably  the  reason  why  that 
procedure  has  proved  successful  here. 
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Figure  1  Effects  of  Inflow  and  Blade  Flapping  on  Vertical  Acceleration 
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Figure  2  Response  to  Collective  Step  Up 


Figure  3  Response  to  Collective  Step  Down  (Case  7) 


Figure  4  Response  to  Collective  Pulse  (Case  9) 
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5.2  Identified  Parameters  (except  rotor  speed) 

The  identified  values  of  all  parameters  not  associated  with  rotor  speed  are  shown  for  all 
cases  analyzed  in  Figures  5  to  10,  while  the  resulting  mean  values  and  standard 
deviations  are  summarized  in  Table  3.  In  the  figures  the  identified  values  as  well  as 
Cramer  -  Rao  bounds  are  shown,  both  normalized  by  the  respective  a  priori  values.  The 
figures  illustrate  that  although  there  is  significant  scatter  in  some  cases,  the  results  are 
mostly  clustered  about  well  identified  means,  i.e.  the  standard  deviations  are  small 
enough  to  provide  meaningful  confidence  intervals.  The  a  priori  estimates  are  also 
shown  in  the  table,  with  a  distinction  made  between  high  and  low  AUW  estimates  where 
appropriate.  Note  that  the  predicted  values  are  derived  using  a  blade  lift  curve  slope  of  6.2 
and  an  'uncorrected'  Mu  value  of  0.85 

A  comparison  of  predicted  estimates  with  values  identified  using  the  current  model 
structure  reveal  a  number  of  interesting  features.  The  vertical  damping  derivative,  -Zw,  is 
significantly  higher  than  the  predicted  value.  This  may  be  attributed  to  the  simple  rotor 
aerodynamics  representation  and/or  may  reflect  the  neglect  of  any  fuselage  aerodynamic 
effects.  For  the  Z e  derivative  the  identified  and  predicted  values  are  within  one  standard 
deviation  of  one  another  for  high  AUW,  but  the  low  AUW  mean  is  approximately  1.3 
standard  deviations  less,  in  absolute  value,  than  the  prediction.  Nevertheless,  this 
represents  a  discrepancy  of  only  11%.  Note  also  that  there  is  some  scope  for  varying  the 
predicted  parameter  values  since  some  of  the  physical  data,  e.g.  blade  inertia,  are  subject 
to  some  degree  of  uncertainty.  For  both  Zw  and  Zq  the  trends  wiht  AUW  appear  to  bo 
consistent  with  that  expected,  but  statistically  meaningful  conclusions  are  not  possible  due 
to  the  relatively  large  standard  deviations. 

The  identified  values  of  the  flapping  derivatives,  pp  and  Pe,  are  generally  within  one 
standard  deviation  of  their  a  priori  estimates.  The  large  scatter  for  the  low  AUW  cases  is 
due  to  high  values  for  the  two  derivatives  obtained  from  Case  7,as  well  as  to  the  small 
number  of  samples. 

The  mean  value  of  the  inflow  derivative,  U©  is  in  good  agreement  with  the  predicted 
estimate,  assuming  Mu  equal  to  0.85.  From  Appendix  A,  t>o  is  seen  to  be  inversely 
proportional  to  M  y  \ ,  so  that  a  'corrected'  value  for  M  n  of  0.54  would  imply  a  t>e  of  over 
2200,  well  above  the  identified  value,  even  accounting  for  the  scatter.  Further,  the 
identified  value  of  t)w  is  most  closely  approximated  by  the  formula  derived  from 
Equation  3  with  Mu  equal  to  0.85.  Thus,  within  the  constraints  of  the  present  model,  the 
dynamic  inflow  model  of  Carpenter  and  Fridovich  [  10]  appears  to  provide  estimates  of  oq 
and  \)w  consistent  with  the  flight  data. 

With  regard  to  the  scatter  of  the  identified  parameters,  it  is  of  interest  to  compare  the 
computed  standard  deviations  with  the  Cramer  -  Rao  lower  bounds.  While  the  latter  vary 
somewhat  from  case  to  case,  the  average  Cramer  -  Rao  bounds  for  the  present  results  are 
typically  3  to  5  times  less  than  the  standard  deviation  which  characterizes  the  scatter.  This 
is  roughly  in  line  with  the  experience  of  the  fixed  wing  community  where  it  is  common 
practice  to  multiply  the  Cramer  -  Rao  bounds  by  a  'fudge  factor'  of  5  to  10  (Ref.  16).  The 
present  results  confirm  the  usefulness  of  the  Cramer  -  Rao  bound  as  an  estimate  of  the 
accuracy  of  identified  helicopter  parameters,  provided  a  suitable  multiplying  factor  is 
used.  Alternatively,  it  provides  some  indication  of  relative  accuracy  of  a  given  parameter 
from  case  to  case.  For  example,  the  relatively  large  Cramer  -  Rao  bounds  for  vq  in 
Cases  1,  and  3  suggest  a  lower  confidence  in  these  points,  whereas,  for  Zw,  the 
relatively  low  Cramer  -  Rao  bound  for  Case  9  adds  credibility  to  the  somewhat  lower 
value  obtained  in  that  run,  the  only  example  of  a  pulse  input. 
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Figure  8  Identified  Results  for  Pq 
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Figure  9  Identified  Results  for  vw 
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Figure  10  Identified  Results  for  V)0 


Table  3  Means  and  Standard  Deviations  of  Identified  Parameters 


Parameter 

A  Priori 
(High  AUW) 

A  Priori 
(Low  AUW) 

Identified 
(High  AUW) 

Identified 
(Low  AUW) 

Identified 

(All  Cases) 

Zw 

-1.03 

-1.16 

-1.57±.22 

-1. 631.39 

-l.59l.26 

Ze 

-465 

-521 

-431±48 

-462145 

-441147 

Pp 

-557 

-557 

-631±64 

-6391251 

-6331135 

Pe 

658 

658 

695+83 

7191248 

703+141 

Dw 

5.8 

5.7 

6.1±1.3 

7.01.92 

6.411.2 

we 

1409 

1409 

1323±268 

12841166 

1310+229 

zn 

-4.6 

-4.7 

14.6±5 

18.016 

15.715 

Pn 

2.1 

2.0 

22.3+23 

38.4122 

27.7123 

WQ 

14 

13 

-108142 

-123176 

-II315I 

«n 

-0.49 

-0.43 

-2.210.4 

-1.510.5 

-1.910.5 

-41 

-38 

-3215 

-251.2 

-2915 

«Q 

.000093 

.000093 

.000097 

.000085 

.000093 

1.000013 

+.000012 

1.000013 

S.3  Rotor  Speed  Parameters 

Table  3  lists  values  of  several  parameters  associated  with  the  rotor  speed.  While  the 
current  model  has  successfully  accounted  for  the  effects  of  rotor  speed  variations  on 
vertical  acceleration  response,  and  has  achieved  a  good  fit  to  the  rotor  speed  record,  it  is 
clear  from  the  identified  parameters  that  the  model  is  not  physically  correct.  Apart  from 
£2 q  (the  inverse  of  the  rotor  inertia),  and  possibly  £2e,  the  other  rotor  equation  parameters 
differ  considerably  from  their  predicted  values.  Similarly,  the  identified  values  of  Zw, 
Pn,  and  uq,  while  subject  to  considerable  scatter  are  generally  much  greater  than  the 
predicted  values,  i.e.  the  predicted  values  underestimate  the  coupling  between  rotor  speed 
and  vertical  response.  Examination  of  the  rotor  speed  record  reveals  a  period  of  around 
1.5  to  2  seconds.  The  estimated  undamped  natural  frequency  of  blade  lagging  motion  is 
about  0.9  Hz,  which,  for  a  damping  ratio  of  0.7,  translates  to  a  damped  natural  frequency 
of  0.64  Hz,  or  a  period  of  1.6  seconds.  The  suggested  implication  is  that  the  blade 
lagging  motion  may  be  an  important  factor  in  coupling  the  rotor  speed  variations  to 
vertical  response.  This  needs  further  investigation,  but  it  is  clear  that  the  extra  parameters 
introduced  by  modelling  the  blade  lagging  dynamics  imposes  an  additional  burden  on  the 
identification  algorithm. 

While  the  coupling  between  rotor  speed  and  vertical  response  needs  further 
investigation,  it  is  not  expected  that  the  parameters  not  associated  with  die  rotor  speed  will 
be  much  affected,  due  to  the  separation  in  frequency  of  the  initial  response,  which  largely 
determines  the  values  of  those  parameters,  and  the  subsequent  periodic  response. 
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Possible  exceptions  are  Zw  and  vw  which  require  time  for  vertical  velocity  to  develop. 
This  has  been  confirmed  by  matching  a  simplified  model,  without  any  rotor  speed  effects, 
with  the  vertical  acceleration  response.  The  match  shown  in  Figure  1(c)  is  typical  of  the 
results.  Significantly,  the  pulse  response  record,  Case  9,  fails  to  converge  with  the 
simplified  model.  This  is  quite  likely  due  to  the  difficulty  in  separating  initial  and  medium 
term  responses  with  this  input.  Neglecting  this  case,  the  mean  parameter  values  for  the 
remaining  cases  are  similar  to  those  given  in  Table  3,  mostly  lying  within  one  standard 
deviation  of  them.  The  main  exception  is  -Z*  which  at  1.93  is  somewhat  greater, 
probably  reflecting  the  poorer  overall  match  of  the  vertical  acceleration,  and  hence  vertical 
velocity.  Significantly  uq  at  1425,  is  very  close  to  the  previous  value  although  t)w,  at  7.6, 
is  slightly  higher,  most  likely  for  the  same  reason  as  for  Zw.  As  expected,  the  Cramer  - 
Rao  bounds  for  the  simpler  model  are  higher,  by  about  50%,  than  those  for  the  model 
including  rotor  speed  effects. 

6.  CONCLUSIONS 

This  paper  has  illustrated  the  use  of  time  domain  Maximum  Likelihood  parameter 
estimation  techniques  in  developing  a  model  of  rotorcraft  vertical  acceleration  response  in 
hover,  capable  of  reproducing  all  the  features  observed  in  the  flight  data  of  a  Sea 
King  Mk  50  helicopter.  A  linearized  vertical  response  model  has  been  derived  from  first 
principles  to  include  both  blade  flapping  and  inflow  dynamics,  as  well  as  rotor  speed 
dynamics.  It  has  been  demonstrated  that  an  excellent  match  of  the  flight  data  can  be 
obtained  with  a  model  including  all  these  features.  A  total  of  nine  flight  cases  has  been 
analyzed  to  identify  model  parameters  for  comparison  with  the  a  priori  predictions.  By 
making  use  of  information  from  an  order  of  magnitude  analysis  and  by  imposing 
constraints  suggested  by  the  a  priori  model,  the  total  number  of  unknown  parameters  has 
been  reduced  to  manageable  proportions.  In  addition,  a  procedure  to  uncouple  the  rotor 
speed  dynamics  from  the  remaining  system  was  developed  to  take  advantage  of  the 
difference  in  frequency  content. 

The  scatter  in  the  results  and  the  problems  with  convergence  of  the  algorithm  were  due 
in  some  part  to  the  large  number  of  system  parameters,  the  limited  number  of  flight 
response  characteristics  recorded,  noise  in  the  recorded  data,  and  the  non-optimal  input 
shapes  used.  Improvements  in  many  of  these  areas  would  be  desirable  in  future  trials, 
especially  in  the  measurement  of  additional  outputs,  such  as  blade  flapping,  and  possibly 
blade  lagging  angle.  This  would  make  possible  the  relaxing  of  some  of  the  constraints 
found  necessary  in  the  current  investigation. 

Despite  the  experimental  scatter,  a  number  of  meaningful  conclusions  can  be  drawn 
from  a  comparison  of  the  identified  parameters  with  predictions.  In  particular  the  inflow 
parameters  are  consistent  with  the  Carpenter  -  Fridovich  dynamic  inflow  model  with  an 
apparent  additional  air  mass  parameter,  M|  i,  of  0.85.  The  vertical  damping  parameter  is 
consistently  higher  that  predicted  possibly  due  to  simplifications  in  the  aerodynamic 
model  used.  Other  identified  parameters,  apart  from  those  associated  with  the  rotor  speed, 
are  in  reasonable  agreement  with  the  predictions.  The  Cramer  -  Rao  bound  appears  to  be  a 
good  indicator  of  relative  accuracy,  and  greater  weighting  could  be  given  to  those 
identified  parameters  with  lower  Cramer  -  Rao  bounds,  all  else  being  equal.  Poor 
agreement  for  the  rotor  speed  parameters  indicates  a  shortcoming  in  the  underlying 
physical  model.  The  results  suggest  that  blade  lagging  dynamics  may  need  to  be  included 
to  explain  the  interaction  between  the  rotor  speed  and  vertical  acceleration  response. 
Further  investigation  is  required  in  this  area. 
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Appendix  A  -  A  Priori  Expressions  for  Model  Parameters 
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Appendix  A  -  A  Priori  Expressions  for  Model  Parameters  (corn) 
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1  Alternative  expressions  from  Equations  (3)  or  (3a)  respectively 
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